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Abstract 

We explain a method for computing the bulk viscosity of strongly coupled thermal plas- 
mas dual to supergravity backgrounds supported by one scalar field. Whereas earlier investi- 
gations required the computation of the leading dissipative term in the dispersion relation for 
sound waves, our method requires only the leading frequency dependence of an appropriate 
Green's function in the low-frequency limit. With a scalar potential chosen to mimic the 
equation of state of QCD, we observe a slight violation of the lower bound on the ratio of 
the bulk and shear viscosities conjectured in pp. 
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1 Introduction and summary 

Recent interest has attached to the possibility of a dramatic rise in the bulk viscosity ( near 
the cross-over temperature T c of QCD where confinement and chiral symmetry breaking set 
in [2], [3]. Bulk viscosity is accessible via AdS/CFT if we break conformal symmetry, for 
example by adding one or more scalars with non-zero profiles in the bulk. Bulk viscosity was 
studied in a holographic setting in [HE] and subsequent work includes [61 [7J [HI [1] . Recently, 
we showed in collaboration with A. Nellore [9] that if a single scalar is coupled to gravity and 
a potential is chosen for it to mimic the QCD equation of state, then the bulk viscosity rises 
significantly near T c , but not as dramatically as expected for pure Yang-Mills theory based 
on the works [21 [10]. The numerical studies in [9] led to the conjecture that the only way 
to get C/s to diverge in the context of a holographic dual in the two-derivative supergravity 
approximation is for the entropy density s to have an extremum as a function of temperature. 
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Figure 1: A cartoon of the relation between shear viscosity and h\2 graviton absorption, and 
between bulk viscosity and absorption of a mixture of the hu graviton and the scalar (j>. 



The purpose of this paper is to provide an exposition of the methods used for calculating 
the bulk viscosity in finite-temperature gravity duals involving a single scalar, as explored 
for example in [HI [12]. Earlier works [H El El [7[ El [1] have focused on calculating the leading 
dissipative term in the dispersion relation for sound waves. We instead appeal directly to 
a Kubo formula. Our computation comes down to understanding the bulk propagation of 
some linear superposition of gravitons and scalar perturbations. The strategy, as illustrated 
in figure [H is to start with an infinite, static, thermal background, linearize the equations of 
motion around it, and solve these linearized equations in some appropriate approximation. 
A large-wavelength, low-frequency approximation is the appropriate one because ( describes 
an effect in hydrodynamics, which is to say the infrared approximation to the dynamics. The 
absorption probability of low-energy quanta is proportional to £. This is entirely analogous 
to the well-known computation of shear viscosity [13] , but more technically involved because 
of the mixing between the scalar and metric perturbations. We can decouple these pertur- 
bations through a judicious choice of gauge, namely the gauge where the value of the scalar 
is used as a radial variable. In this gauge, we need only solve a single ordinary differential 
equation, given in equation (pal) . The bulk viscosity to shear viscosity ratio can then be 
extracted from the solution of f)9ap by using equation (j4ip . which is our main result. In 
examples, a conjectured bound [1] on the ratio £/r) of bulk to shear viscosity is observed to 
be satisfied in most but not all circumstances. 

The structure of the rest of this paper is as follows. In section [2] we explain our method 
for computing the bulk viscosity from AdS/CFT. Section I2TT1 summarizes the relevant grav- 
itational backgrounds. Section 12.21 shows how to perturb the backgrounds in a rotationally 
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invariant fashion. Section 12.31 reviews the treatment of metric perturbations that control 
the shear viscosity, explaining in particular how to extract the imaginary part of the re- 
tarded Green's function of an off-diagonal component T 12 of the gauge theory stress tensor. 
Section 12.41 adapts this treatment to the case of bulk viscosity, using equations derived in 
section 12 .21 Section 12.51 detours once again to a review of the shear viscosity computation, 
showing how to solve the linearized equation of motion in the small uj limit and giving the 
final prescription for extracting 77. Section 12.61 explains the analogous computation for the 
bulk viscosity. In section [3] we compute the bulk viscosity for specific backgrounds. For the 
Chamblin-Reall backgrounds we obtain an analytic result for (/r] in section IBTTj where we 
also give a different derivation of the same result using a Kaluza-Klein reduction argument. 
In section 13.21 we show some results of numerical computations of the bulk viscosity in a 
particular class of backgrounds whose equation of state is similar to that of QCD. For com- 
pleteness, we provide derivations of the Kubo formulas relevant to our results in appendix [XI 
In appendix [B] we clear up a subtlety in the computation of two-point functions related to 
our choice of gauge, and in appendix [C] we justify the prescription for the computation of 
the retarded two-point function that we used in section [2j 

2 Computation of viscosity using AdS/CFT 
2.1 The background geometry 

We consider backgrounds with a scalar coupled to gravity. The relevant action is 





and the corresponding equations of motion are 



n<t> = v'(<j>) 




= T, 



(2) 



where the stress tensor for the scalar is 



T, 




(3) 
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The backgrounds of interest have the form 

ds * = e 2AW [- h ( r )dt* + dx 2 } + e 2B{r) ^- $ = cj>(r) . (4) 

h[r) 

The choice of radial variable r is arbitrary: reparameterizing it leads only to a different 
choice of B. A convenient choice for the backgrounds we will study is r = $, namely the 
unperturbed value of <f)(r). Plugging (jlj) into ([2]) leads, with this choice, to 

A" - A'B' + - = 
6 



h" + (4A' - = 
6A'h' + h{2AA' 2 - 1) + 2e 2B V = 

AA' - B' +-- - —V = 
h h 



(5) 



where primes denote d/d<&. Using a numbering scheme (t, oc . ■ oc ■ it} — (cc . 3j . OC ■ OC ■ X 
the first two of these equations come from the 00 and 11 Einstein equations; the third comes 
from the 55 Einstein equation; and the last comes from the scalar equation of motion. There 
is typically some redundancy in equations obtained from classical gravity, with or without 
matter. In the case of (jSJ), the redundancy is that the $ derivative of the third equation 
follows algebraically from the four equations listed. 



2.2 Rotationally invariant deformations 

In order to extract transport coefficients, we must calculate the u — > limit of two-point 
functions of the schematic form 

G(u) = J dtd 3 xe luJt (O(t,x)O(0,0)) . (6) 

As we explain in more detail in appendix |A] when calculating the bulk viscosity, the operator 
O may be taken to be = Xn + T 22 + T 33 . Because this combination preserves the 50(3) 
symmetry of spatial rotations, and because we integrate in over spatial separations x in 
an 50 (3)-symmetric fashion, we should be capable of framing the entire calculation of bulk 
viscosity in an 50 (3)-invariant fashion. Imposing this symmetry on the metric ansatz forces 
all components of g^ v to depend only on t and $, and it also forces gu = (722 = 533 and the 
vanishing of g i0 , g i5 , and g^ for i ^ j. Gauge freedom allows one to stipulate that $ doesn't 
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change at all, and also that gos = 0|_| To summarize, the perturbed metric takes the form 
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where 



9oo 



9n 



955 



i 2A h 
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H, 



55 



A, B, and h are functions of $ only which satisfy the equations (jSJ); H 00 , Hn, and H 55 
depend on t and <3>; and A is a formal expansion parameter. The reason for the precise form 
(JBJ) is that it leads to a relatively simple expression for s/—g. 

In the presence of the ansatz the non-trivial components of Einstein's equations are 
for [iv = 00, 11, 55, and 15. The scalar equation of motion is also a non-trivial constraint on 
the metric components. These five equations may be formulated in terms of <7ooj gn, and #55, 
but their explicit form is complicated and would not add clarity to the current exposition. 
Let us instead pass directly to the linearized approximation, where we keep terms only up to 
O(A) in the equations of motion. Let us also assume harmonic time dependence by replacing 



-iujt 



H^ u (<&). Because of some redundancy in the five equations of motion, 



similar to the redundancy of (jSJ), one winds up with only three independent equations, 



1 Stipulating that $ doesn't change and 1705 = is a convenient alternative to the more standard axial 
gauge, where $ is allowed to change but gos = and 355 is unperturbed. The reason that two perturbations 



may be set equal to zero is that S'0(3)-invariant diffcomorphisms 



£ M are specified by two 



independent functions of t and $, namely £° and £ 5 . There is a subtlety related to our choice of gauge that 
we clarify in appendix [Bl 
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(9c) 



These equations are usually too complicated to solve in closed form; indeed, the same is 
true of the equations (jSJ) that determine A, B, and h. However, a combination of series 
expansions near the horizon (where h has a simple zero as a function of $) and numerics can 
be used to find fairly precise solutions. The terms proportional to u 2 may be ignored except 
very close to the horizon. The horizon boundary conditions appropriate for computing a 
retarded Green's function are that only an infalling wave is allowed at the horizon. 

With a solution to (I9*aj) - (l9"c]) in hand that satisfies appropriate boundary conditions at 
the horizon, one still has to evaluate the on-shell action in order to find Green's functions, 
using the prescription of [HI [To*] . In fact, there are two subtleties involved in this evaluation: 

• The Green's functions of interest are complex, whereas the action is real. The transport 
coefficients are read off from the imaginary part of the Green's functions, which is 
closely related to a conserved Noether current in a complexification of the supergravity 
action. 

• Boundary terms must be added to the action (pQ) in order to obtain a Lagrangian which 
depends only on the perturbations and their first derivatives [16J. Additional 
boundary terms are often needed to cancel divergent parts of the resulting Green's 
functions. 

2.3 Shear perturbations 

In order to illustrate and resolve the two subtleties mentioned at the end of the previous 
section, it helps to start by rehearsing the now-standard case of the shear viscosity. The 
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metric, including perturbations, needed in order to study the shear viscosity is 
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(10) 



where, as usual, Hi 2 is a function only of t and and we have persisted in using the 
gauge r = $. This should be compared to ©-([8]) in the bulk viscosity case. At the level of 
linearized perturbations of the equations of motion, or quadratic perturbations of the action, 
H12 decouples from all other metric perturbations. It is straightforward to show, again to 
linear order in A, that the only non-trivial equation of motion in the presence of the ansatz 
f fTUj) is the 12 Einstein equation, which reads 



H'[ 2 + ( 4A' - B' + 
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where we have again substituted H 12 (t, $) 



-iuit 



iJ 12 ($). The result f[TT|) in the shear 



viscosity computation is parallel to fl9a|) -fl9cl) in the bulk viscosity computation. Assuming 
that we have found a solution (at least an approximate one) to f fTTj) . we must now confront 
the evaluation of the on-shell action to 0(X 2 ). The first difficulty is that the substitution 
Hi 2 (t, $) — * e~ luJt Hi 2 (&) doesn't make sense outside the context of a linearized equation. 
So let's instead keep H\ 2 as a general (real) function of t and $. The next difficulty is 
appropriately fixing total derivative terms. Starting from without any improvements or 
boundary terms added, one finds 
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16ttG 5 
£ = t + d t t + 
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£ = 

2h 



rt _ Z J2A+B TT TT 

L — — — e H 12 n 12 
h 

= 2he iA ~ B H l2 H[ 2 + hA'e 4A ~ B H 



(12) 
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We will refer to C as the "unimproved" lagrangian and to £ as the "improved" lagrangian. 
The improvement terms are total derivatives dtC 1 and <9$£*, and they come in part from the 
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well-known Gibbons-Hawking term involving the trace of the extrinsic curvature integrated 
over the boundary of spacetime. The form of the improvement terms is largely fixed by 
demanding that £ should depend only on H\ 2 and its first derivatives. But this requirement 
doesn't constrain the second term in as shown in ffl2l) : indeed, any multiple of Hf 2 could 
be added to and correspondingly one would wind up with terms proportional to Hi 2 H' l2 
and H\ 2 in the improved lagrangian. The particular choice of we made in (fT2"j) obviously 
leads to the simplest form of C, but in the treatment of the bulk viscosity it is not so obvious 
how to choose Let us therefore consider a slightly generalized "improved" action and 
lagrangian: 

S = / d 5 x L 

167rG 5 J ^ 

2h 12 2 12 2 12 12 

£ evidently differs from the improved lagrangian in (Tl2l by the total derivative <9$ (|Gi? 2 2 ). 
Here G is an arbitrary, smooth, real function of $. The idea (modulo the issue of real versus 
complex perturbations) is that evaluating S by plugging in a solution to the linearized 
equation of motion ffTTT) gives a two-point function of the operator dual to if 12, namely X12. 
But as long as G is not fixed by some first-principles consideration, this Green's function will 
be ambiguous! Fortunately, we will soon see that this ambiguity affects only the real part of 
the Green's function, and we want the imaginary part. This is not to say that G cannot be 
determined: it is closely related to counterterms in the prescription for computing Green's 
functions which are widely studied: see for example [171 US EE] • The point is that we do not 
need to determine it. 

Having side-stepped the more subtle aspects of boundary terms, let us move on to the 
issue of obtaining a complex-valued Green's function from a real-valued action. It helps to 
pass to the following real- valued lagrangian for a complex- valued field h 12 = h 12 (&): 

C c = ^-e 2A+B \h 12 \ 2 - he 4A - B \h' 12 \ 2 + G'\h l2 \ 2 + G(h 12 h? 2 + h\ 2 h' l2 ) . (14) 

£c is constructed so that the equation for /i 12 obtained from it is identical to the one 
obtained for H12 from the improved lagrangian in (Tl3l) with the assumption of a harmonic 
time dependence e~ lu)t — that is, it is the 12 Einstein equation, fill I) . One may re-express Cc 
as a total derivative plus terms that vanish upon use of the equations of motion: 

Cc=d * J+h Awi~^dhg) (15) 
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where 

J = -he 4A ~ B h* 12 ti 12 + G\h 12 \ 2 . (16) 

Up to a real normalization factor and the ambiguity in G(<&), the limit of J as one approaches 
the conformal boundary is the quantity to be identified as the retarded Green's function. 
This claim essentially follows the discussion of [20J, which was justified using a first-principles 
application of the prescription of [TH [T5] in a Schwinger-Keldysh formalism in [2T]. In 
appendix [Cj we provide a different derivation of this result. Because we only care about the 
imaginary part of the Green's function for purposes of computing transport coefficients, we 
define 

T = - Im J = -he 4A - B {h* 12 h' 12 - hX 2 h 12 ) . (17) 

JF coincides precisely with the conserved Noether charge associated with the U(l) symmetry 
of ([Hj) under a phase rotation h±2 — > e ld h\2. It has the physical interpretation of the number 
flux of gravitons in the radial direction^ The sign of JF was chosen so that JF > corresponds 
to particles falling into the black hole. T is evidently the Wronskian of h\i and h\ 2 , multiplied 
by the factor required by Abel's identity to make it constant when the equation of motion 
( ITT]) is obeyed. As such, its form up to an overall multiplicative constant could have been 
guessed directly from the equation of motion (ITT]) without detailed considerations regarding 
the on-shell action. 



2.4 The number flux for rotationally invariant perturbations 

In analogy to (fTTj) . it is natural to guess that the imaginary part of the retarded Green's 
function for is proportional to a flux proportional to the Wronskian of a complexified 
solution tin with its conjugate: 

The $- dependence of the prefactor can be deduced from Abel's identity applied to (JSaj) , but 
in our current approach, the overall constant factors in (fT8|) must be fixed by a calculation 
of the on-shell action. We describe such a calculation in the next paragraph. The algebra 
gets a little complicated because there are three fields involved, namely H 00 , if n , and H 55 . 
However, every step is analogous to the foregoing discussion of the on-shell action in the 
shear viscosity calculation. 

2 Indced, calculations of imaginary parts of Green's functions using conserved number fluxes were com- 
monplace before the prescription [TU [15] was suggested: see for example J22J [23l EH [25] . 
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The action §B) in the presence of the ansatz Q takes the form 
1 
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(19) 



(20a) 



(20b) 



(20c) 



(20d) 



and G is a symmetric matrix which may depend on $ in an arbitrary but smooth way. In 
principle, G should be determined by insisting that the full action has only first derivative 
terms; however, as in the shear viscosity case, we do not require the explicit form of G. 
Likewise, we do not need to know the explicit form of the improvement terms L l and 
The real-valued lagrangian analogous to ([141) for complex-valued fields h = h(&) is 



C c = h*' T mh' + h* T kh + h*' T bh + h* T b* T ti 



(21) 



where 



m = M 



k = u 2 M tt + M + G' b = M* + G. 



(22) 



10 



Just as in (1T51) . we may express Lq as a total derivative plus terms that vanish on-shell: 



Cc = d^J + h 




dh* d<& dh*' 



) 



(23) 



where 



J = h* T (mh' + hh) . 
The conserved Noether charge associated with the symmetry h — > e %e h is 



(24) 



T = - Im J = - \h* T (mh' + hh) - (h*' T + h* T h T )h 



(25) 



The form (125]) appears to be distinct from (|T8l) . But one can see from (pal) -(|9c |) that all 
quantities appearing in T as defined in (1251) can be eliminated in favor of h' lx and hu. It is 
straightforward to check that the result of such an elimination is precisely the form quoted 
in f|T8l) . As in the case of the shear viscosity, T has the physical interpretation of a conserved 
number flux of particles falling into the horizon. The particles in question are a mixture of 
metric and scalar perturbations, though as we have seen, a gauge choice is possible that sets 
the scalar perturbations to zero. 

2.5 Low-frequency limit of shear perturbations 

With the conserved number flux in hand, we should be able to extract the imaginary part 
of the retarded Green's function given a solution of the equations of motion. As before, it is 
illustrative to start by reviewing the shear viscosity computation. The linearized equation 
of motion ffTT]) can be solved for uj = 0: 



where a\2 and byi are integration constants and we have assumed that $ — >• at the conformal 
boundary. Because h has a simple zero at the horizon, $ = the second term diverges 
logarithmically there. In the strict limit uj = 0, this term is disallowed. For very small uj, we 
may use a matching procedure which is analogous to boundary layer theory and well known 
in the literature on absorption by black holes: see for example [26]. The trick is to find the 
leading behavior of solutions for non-zero uj just outside the horizon: 



uj = solution: 




(26) 



$h solution: 



+ q 2 (^-<i>)^ /47rT , 



(27) 
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where the Hawking temperature is 



-h'($ H ) e A{$ H )-B{$ H ) ^ (28) 



47T 



and we have assumed that > and h — > 1 as $ — > 0, i.e. at the conformal boundary. 
Standard horizon boundary conditions, associated with the enforcement of causality, are to 
suppress the solution corresponding to gravitons coming out of the horizon: that is, c\ 2 = 0. 
The matching procedure is to expand fl26|) around the horizon and ( 128]) around uj = and 
extract relations among a 12, 612, and c{ 2 by comparing the expansions. Briefly, the result is 



e ~4A(^ H )+B(^ H ) 



hi2 ~ a 12 + 612 Z77^ \ lo g( $ H - $) ~ c i2 



1 — 10g($f7 - $) 



(29) 



From (1291) we find, to leading order in small uj, 

= c{ 2 b 12 = tue 3A ^c^ 2 . (30) 
The flux is most easily computed in the limit $ — > from (I27|) : 

^ = we 3A(* H )| c -|2 ; (31) 

still at leading order in small uj. 

The imaginary part of the retarded Green's function is 

ImG * M = -l6^ (32) 

where 

G r {uj) = -1 [ dtd 3 xe l ^e(t)([T 12 (t,x),T 12 (0,0)]) , (33) 



and T is computed assuming a i2 = 1. The shear viscosity is 

1 P 3A(# H ) 

77 = - lim - Im Gji(w) = — . (34) 

c<J 107rG5 

One may confirm the relation 7]/s = l/4ir at this point simply by noting that the entropy 
density is given by s = e 3A ^ H ^ /AG^. It is also worth pointing out that we never needed the 
explicit result for fe 12 quoted in (130|) . because we computed the flux T at the horizon, not at 
the conformal boundary. In fact, we only needed the uj = solution which is regular at the 
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horizon and approaches one at the boundary; c 12 is then given exactly by the value of Hi 2 
at the horizon. 



2.6 Low frequency limit of rotationally invariant perturbations 

The imaginary part of the retarded correlator of T { % can be obtained in a similar manner 
by solving (I9al) . The u = equation doesn't seem to have an analytic solution in general, 
so it is usually necessary to proceed numerically. The boundary conditions that we need to 
impose are hu(0) = 1 and the fact that hu should be regular at $ = As in the shear 
viscosity case, this solution needs to be matched onto the boundary layer at $ « 
Close to the horizon, the two solutions of fl9al) behave as 

$ « $h solution: h n = c+($ H - $)^/ 47rT + c u (<t> H - $)~^/ 47rT , (35) 



where T is the Hawking temperature as in (J28J). Infalling boundary conditions mean setting 
c u = 0. For small a;, 

r/ti 

(36) 



/in ~ c n 



log($fl- - $) 



4ttT 



so to leading order in we have hu « c^. We can therefore find by just evaluating the 
a; = solution at the horizon. 

With the value of c n in hand, we can find T by plugging ( 1351) into ( 1181) and keeping 
only the leading order term in — $. The subleading terms do not give corrections to 
J 7 , because T is conserved so it should be independent of $. The fact that upon plugging 
(1351) into (TIB]) we also get terms suppressed by various powers of — $ just shows that the 
asymptotic form (135j) is not exact and should in fact be corrected in such a way that the 
resulting T stays independent of $. We obtain 

r =^ mH) M»f- (37) 



Given J 7 , we can extract the imaginary part of the retarded Green's function of T { 



throug 



ImGRM = -i6iv (38) 



3 The coupling of the metric to the stress tensor is ^H mn T mn . So when we set Hu = H 22 = H33, the 
dual operator is . But when we set H\<x = the operator we couple to is T\ 2 . That is why there are 
explicit factors of i in ([55)) but not in 
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where 

G R {u) = -i J c/trf 3 xe^^(t)([ir/(t,f),iT fc fc (0,0)]). (39) 
The bulk viscosity is then given by the Kubo formula (see appendix [A]) 

C = ~\ lim-ImG^H = J |c n | 2 , (40) 

where we have used s = e 3A ^ H ^ / AG^. Since Einstein's equations §5$) for the background 
imply A'{$ H ) = -V($h)/3V'($ h ), and rj/s = l/4vr, equation (@DD yields 



,2 



c n l • (41) 



:2 1^11 



In [TT] it is noted that the speed of sound in the background PJ can be approximately 
expressed as 

i _ i vw (42) 

Combining this with (|4"IT) we obtain 



„ 2 U- C -J |C "' - (43) 

It is interesting to compare this expression with the conjectured bound on (,/r) of [I] stating 
that all 3 + 1-dimensional finite-temperature field theories should satisfy 



^Hs-^J' (44) 

We will refer to this inequality as the viscosity ratio bound. Using (l4"3l) we note that ( 1441) is 
roughly equivalent to (c^j 2 > 1 — roughly because (143]) itself is an approximate result. 



3 Examples 

3.1 Chamblin-Reall backgrounds 

The simplest application of the result (1411) is to Chamblin-Reall backgrounds [27], which are 
gravity solutions following from the action ([1]) with potential 

V(<j>) = V e*, (45) 
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where Vq < and < 7 < \/2/3 are constants. In a gauge where r = $, the coefficients of 
the Chamblin-Reall metric take the form 

24 ( 2$\ 2B 8 - 37 2 , _ 

e lA = exp ( -— ) e 2 = — ^ exp(-7$) 

(46) 



h = 1 — exp 



37 / 67 2 Vc 
-3 7 



6^ 



An interesting feature of the Chamblin-Reall solutions is that the entropy density has a 
power-like dependence on temperature, 

s oc T^ 2 " 3 ^ . (47) 

If this background were dual to a field theory, equation (1471) would imply that the speed of 
sound in such a theory is constant and given by 



c 



2 



rflogT _ 1 
d log s 3 



2 



(48) 



The first equality above is true in any field theory at zero chemical potential. 

It is straightforward to check that for a Chamblin-Reall background the coefficient of hu 
in (15al) vanishes when uj — 0. Consequently, hu = const is a solution to the u = equation. 
Since the other linearly independent solution diverges logarithmically at the horizon, and 
since one requires hu — > 1 in the ultraviolet, one concludes that c{ x = hu(<$>H) = 1 in this 
case. The ratio of bulk to shear viscosity can then be easily found from (|4ip to be 

£ -.2 n r „2 



, ^H§- c -J- (49) 

This shows that, as was already pointed out in p], the Chamblin-Reall backgrounds saturate 
the viscosity ratio bound (]44j) . 

There is a more conceptual way of understanding (I49p that comes from the fact that 
the Chamblin-Reall metric (]46p can be obtained through a Kaluza-Klein reduction from a 
.D + 1-dimensional anti-de Sitter space AdSo+i- (Related arguments are presented in [El II]-) 
The Kaluza-Klein reduction is explained in [H], but the details of this reduction will not be 
relevant in the following argument. All we will need is the relation between 7 and D: 

D= S -^. (50) 
2 - 37 2 V ; 
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In any field theory in d spacetime dimensions, the bulk and shear viscosities are related to 
the two-point functions of the stress-energy tensor by the generalization of the Kubo formula 
§8}: 

V [SikSji + 6«5 jk - -T^-rSiAi ) + (SijSu = lim - / rf 3 xrfte^^(t)([T i ,(t,x),T fci (0,0)]) , 
\ a — 1 J u— >o to J 

(51) 

where all indices denote spatial components. Setting i — j and k = I and summing over i 
and k running from 1 to 3, we obtain 

3 

9Ccr = lim - / d 3 xdte wt 9(t) V ([T«(t, x), T kk {0, 0)]) CR (52) 

J i,fc=l 



in the Chamblin-Reall case and 

?/cft ( 6 - -r^-r^j = lim - / d 3 x e^*0(t) V f), T fcfc (0, 0)]) C ft , (53) 

in the CFT/) dual to the anti-de Sitter space AdS^+i from which the Chamblin-Reall solution 
was obtained. In (153]) we used the fact that Ccft = by scale invariance, so there is no bulk 
viscosity contribution. 

Because field theory two-point functions can be computed from taking functional deriva- 
tives of the on-shell gravitational action, the two-point function in the right-hand side of 
( 152]) differs from the one in (1531) by the ratio Gd+i/G^ of the Newton constants in D + 1 and 
in 5 dimensions. Consequently, 

E?,fc=i([ T »(^ x),T kk (0, 0)])cr _ G D+1 _ s C r Vcr (t .,, 

E', fc =l([^u(^^)^fcfc(0,0)])cFT " G b " S C FT ~ VCFT ' 

In the last equality we used the fact that rj/s = l/An in both theories. Combining f[5"4"]) with 
([52]) and ([53]), it follows that 

CCR 2 (\ - -L-) . (55) 



?7cr V 3 D-l 

Since (1481) and ([50]) imply that in the Chamblin-Reall background (? s = 1/(D — 1), one can 
immediately recover ([413]) . If we now set D = 6 we obtain (/rj = A/15. This result was 
obtained in [7] through different methods. 
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Figure 2: s/T 3 as a function of T/T c for two potentials of the form (|56|) with 
{7 0.606, b ~ 2.06, A ~ 3.93} and {7 « 0.606,6 « 1.503, A « 3.61}. For each curve, 
T c is defined to be the ordinate of its inflection point. These potentials were obtained in 
[TT| [9] where their equations of state are shown to mimic that of QCD. 

3.2 Numerical results on QCD-like models 

With (14 ip in hand, we can compute (/s numerically for any choice of the scalar potential 
V{4>) over a range of temperatures. In this section we will present results for a specific form 
of V((j)) chosen so that the equation of state of the black hole resembles that of QCD [TT| [9] : 



Following the nomenclature of [9] , we describe black holes derived by coupling a scalar with 
this potential to gravity as Type I black holes. The quantity L is the radius of curvature of the 
asymptotic AdS geometry. We choose to work in units where L — 1. We choose 7 ~ 0.606 in 
(]56p in order for the squared speed of sound in the IR to be approximately (? s ~ 0.15, a value 
broadly consistent with hadron gas phenomenology. Similarly, the parameter b appearing in 
(|56p can be adjusted so that the dimension A of the field theory operator dual to the bulk 
field matches that of tr in QCD at a particular scale. As explained in [9], A w 3.93 at 
an energy scale Q = 3 GeV, which gives b ~ 2.057. If the matching is done at Q — 1 GeV, 
one obtains A ps 3.61 and b ~ 1.503. Figures [2] and [3] show the equation of state and £/s 
corresponding to (1561) for these two choices of b. Figure [3] should be compared to figure 1 
of [3], representing the ratio £/s for QCD as computed from lattice data. Similarly to what 



12 
1? 



cosh (70) + bcj) 2 . 



(56) 
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Type I BH, A=3.93 



Type I BH, A=3.61 




Figure 3: The ratio (/s (dotted) for two potentials of the form (1561) . whose equations of 
state were plotted in figure [2j The solid curves show the Buchel bound (jHJ), and a violation 
is evident when T is below T c for both potentials. 

was found in [3] , (/s does increase as T approaches T c from above. However, this quantity 
rises only up to C/ S l m ax ~ 0-06 around T c , in contrast with the sharp rise exhibited in [3]. 

It is worth mentioning that the values of (/s obtained in these examples violate the 
viscosity ratio bound (jHJ) for a range of temperatures below T c . As can be seen from 
figure [3j in the case A pa 3.93 the bound is violated for 0.541 < T/T c < 0.900, while in 
the case A pa 3.61 it is violated for 0.345 < T/T c < 0.737. Even though the violations 
are smaller than 5% for the potential with A pa 3.93 and smaller than 3% for the one with 
A pa 3.61, they are well within the margin of numerical error. But although the potential 
fl56|) is similar to those encountered in gauged supergravity, its precise form was made up. So 
the possibility remains that the bound f|44|) could be true in all potentials actually derivable 
from string theory. 

By solving f[9"al) for finite values of uj and computing the flux JF as given in (ITgj) . we can 
also evaluate the spectral density function given by 

p(uj) = --ImG R {u) = -^—. (57) 

7T 167T Z G5 

As is explained in [3], in QCD p ~ a> 4 for large values of u. In our model, preliminary numer- 
ical studies at large u suggest a power-law behavior ImGfl(w) oc uo u but with v depending on 
temperature. For instance, for the potential of the form (156]) with {7 pa 0.606, b = 2.06, A pa 
3.93} we found that v pa 2.8 for T/T c pa 1.05 and that v pa 3.2 for T/T c pa 1.45. These values 
of v were obtained from power-law fits in the range 5 < uj/T < 20, which may be too small 
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a window to accurately capture the high-u; behavior. 
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A Kubo formulas for bulk and shear viscosity 



For completeness, we now outline the proof of the general Kubo formula for shear and bulk 
viscosity 



V ( S ik 5ji + 5u5 jk - -SijSki) + CSijSu = - lim — lmGf iM (uj) 

3 ) (58) 

Gf jM {^) = -i / rf 3 x^e^(t)([T y (t,f),T fci (0,0)]), 



from which one can obtain ( 1341) by setting % = k = 1 and j = I = 2, and equation ( 140|) 
by tracing over ij and fci. All indices in ( 158|) run over spatial directions only. The proof 
consists of two parts. In section PV. II we apply Kubo's linear response theory to express the 
stress-energy response to metric perturbations in terms of the retarded two-point functions 
of Tij. In section IA.2I we extract the shear and bulk viscosities by comparing the results 
in section IA.1I to linearized hydrodynamics. The results included in this section are not 
new: Kubo's linear response theory was developed about fifty years ago in [281 129], and the 
stress-energy response to metric fluctuations in linearized hydrodynamics is explained, for 
example, in [30 J. 

A.l Kubo's linear response theory 

We start with a finite-temperature quantum mechanical system whose Hamiltonian is 

H = H + H' H'= -\05(t) , (59) 

where O is some time-independent Schrodinger operator and A is a small parameter. In 
units where H — 1, the equation of motion for the density matrix, 

«'^ = [*M. ( 6 °) 

can be solved through time-dependent perturbation theory: 

p(t) = e- lHot p e iHot + Sp(t) 5p(t) = ie~ iHat {\(D, p ]e iHot 9(t) (61) 

where po is the unperturbed density matrix. Note that the above equations are all written 
in the Schrodinger picture where p evolves in time, but O doesn't. 

The expectation value (V(t)) of a different operator V at a later time t can then be 
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computed from 



(V(t)) = tr {Vp(t)} = tr {Pe- iHot pe lHot } + 6{V(t)) 

(62) 

6(V(t)) = iti{e- iHot [\0,p ]e iHot P} . 
Using the cyclic property of the trace and passing to the Heisenberg picture, one obtains 

5(V(t)) = i6(t)tr{p[V(t),\O(0)}} = i6(t){[V(t),\O(0)}) . (63) 

The factor of 8(t) in the above formula can be understood from the fact that the perturbation, 
which occurs at t — 0, can only affect the expectation value of V at later times. Similarly, the 
commutator between V(t) and O(0) arises because it is a measure of the degree to which the 
perturbation can affect V(t). In particular, if V(t) and O(0) commute, then the expectation 
value of V at time t is unaffected by the perturbation, so (5V(t)) = 0. 

The result (1631) can be applied to a field theory on R 3 ' 1 whose lagrangian density is given 

by 

C = C + hi kl T kl , (64) 

where £ is the unperturbed lagrangian density, hki denotes perturbations of the metric, 
and T kl is the stress-energy tensor. Again, indices k and / run only over spatial directions. 
Taking hki(t,x) = hkiS 3 (x)S(t), and using fl63l) with XO = hki5 3 (x) and V(t) = Tij(t,x), one 
obtains 

SiT^x)) = U(t)(\T l3 (t,x),h kl T kl (t,0)}) ■ (65) 

Because small perturbations add linearly, we can extend fl65|) from perturbations localized 
at x = and t — to spatially homogeneous perturbations where hij(t, x) = hij(t): 

5(T, l3 (t,x)) = l -J dt'd 3 x' 9(t - t')([T, tl (t,x),T kl (t',x')])h kl (t') . (66) 

In Fourier space, this equation reads 

SiTyfag)) = ~h u {u)Gf l (u) , (67) 

where Gf jkl (ui) was defined in (158|) . 
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A. 2 Spatially homogeneous perturbations in linearized hydrody- 
namics 



In a curved background, the stress-energy tensor of a viscous fluid with velocity field 
energy density e, pressure p, and bulk and shear viscosities ( and r\ is given by 



T» v = (e + p)«V + pg^ - P^P"? 



(68) 



where P^ u = + is a projection operator and V M denotes covariant differentiation. 
We are only interested in perturbing a fluid at rest on R 3 by changing the space-space 
components of the metric to gij = rjij + hij. In this case, the velocity field = (1,0,0,0) 
doesn't change to first order in because of parity symmetry. It can be checked that to 
lowest non-vanishing order in the perturbation, the change in the space-space components 
of the stress-energy tensor is given by 



1 

3' 



1 

2 ( 



(69) 



where K = —Vdp/dV is the bulk modulus. In frequency space, this relation reads 



pSfS/ - l -K5 l3 5 kl 



-hkiiuj) 



V [ 5*5/ + 6*6? - ) + CM 



ki 



ckl 



(70) 

Because linearized hydrodynamics is expected to be valid in the limit of slowly-varying 
perturbations, one should compare this formula to equation (jBTj) . Equation (158!) follows 
immediately. 



B Gauge transformations and dual operators 

It is standard in AdS/CFT to assert that the leading behavior of the H mn component of 
the metric perturbations translates into a deformation of the field theory lagrangian by the 
operator ^T mn , provided that one works in "axial gauge" where H 5m = 0. But we work in a 
gauge where 5(f) = 0, and as a result H 55 is not 0. The goal of this section is to justify the 
fact that, despite this gauge choice, the operator dual to H u = H 2 2 = H 33 is still given by 
2^ 1 ) as would be the case in axial gauge. Most of the discussion will rely on the following 
asymptotic behaviors of various metric components near the boundary of AdS (see |llj). 
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written in a gauge where r = $: 

A « — - — log $ B « - log $ /i«l 

A-4 (71) 

50 « a$ # mn « i? mn . 

Here _R mn and a are constants, and A is the dimension of the operator 0$ dual to 0. 

Under an infinitesimal gauge transformation, the metric and scalar perturbations trans- 
form as follows: 

$9mn > &9mn &9mn ~t~ ^ ' m£,n ~t~ ^ ' n^Lm (72) 
50 _ 50 = 50 + f m m $ , 

where £ m = £ 1 , £ 2 , £ 3 , £ 5 ) is the vector field parameterizing the gauge transformation. In 
our coordinates r = $, so 50 = 50 + £ 5 . Let's assume that initially 50 = and 

6g mn = diag {-he 2A H 00 , e 2A H n , e 2A H lu e 2A H lu e 2B H 55 /h} . (73) 

We would like to perform a gauge transformation to axial gauge, where 

Sg mn = diag {-he 2A H 00 , e 2A H n , e 2A H u , e 2A H lu o} , (74) 

and where 50 = £ 5 does not necessarily vanish. The symmetries of the problem imply that 
= £ 2 = £ 3 = 0, as well as £° = £°(t, $) and £ 5 = £ 5 (t, $). To achieve the above gauge 
transformation, we need to solve the following two equations, coming from the 55 and 50 
equations in (172]) . respectively: 

e 2B-2A C^) 

It follows from (ITTj) that near the boundary (i.e. at small $), 5' = —1/$ becomes much 
larger than h'/2h, so the first equation in f!75|) has the asymptotic solution 

50 = £ 5 ~ C(t)$ + non-homogeneous solution , (76) 

where C(t) is an integration constant that may depend on t. By adjusting this constant, one 
can therefore construct a gauge transformation between ( 173]) and ( |74l) for whichE S(f> = o($) 



4 We employ here the "little o" Landau symbol: f(x) = o(g(x)) as x —* iff lim^o f(x)/g(x) = 0. 
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as $ — > 0. In addition, the 00 and 11 equations in (1721) give 

H 0O = H O0 + 2A'6^+^-6i + 2d t ^ / x 

h (77) 

H n = H 11 + 2A'5(f), 

which, together with (1751) and (17T1) . also imply that H 00 = H 00 + o($) and i7 n = if n + o($) 
as $ — > 0, provided that A > 2. From 5(f> = o(<3>), together with H Q0 = H Q0 + o($) and 
Hu = Hu + o($), it follows that the operator dual to Hn is the same in a gauge where 
5(f> = and in axial gauge: it is given by \T i l . 

Note that we might equally have chosen a gauge where 5(f> ~ a<3> as $ — ► with non-zero 
a. In this gauge, H 00 = H 0Q + + °( < ^ ) ) an d ^Tli = -Hn + A^i a + °( < ^ ) )- The contribution 
to the deformation of the lagrangian proportional to a is then 

^oc^xirz + ev (78) 

This expression vanishes upon using the fact that = —ftO^ = (4 — A)0^. 



C Prescription for computing the retarded correlator 

The prescription that we used for computing retarded correlators (equations ( l32l and ( |38f) ) 
was proposed in [20], and it was justified in [2T] by using Schwinger-Keldysh contours. In 
this section we will show how this prescription can also be justified through the fact that 
it provides an analytic continuation of the corresponding Euclidean correlator computed 
according to the prescription of [T4"l [T5] . This fact was checked in [20], but a more general 
argument can be made. 

For a bosonic operator O, the Euclidean correlator is defined by 

G E (u n , k) =J P drJ d 3 x e -*#*-<^) (T e O(t, 2)O(0)) , (79) 

where T E denotes Euclidean time-ordering, r is the Euclidean time, and u n = 2irn/f3 are 
the Matsubara frequencies. The Euclidean correlator is defined only at the Matsubara fre- 
quencies u n because in position space the time-ordered two-point function appearing in (1791) 
is periodic with period (3 = 1/T. G E (u n , k) is related to the retarded correlator Gn(uj,k) 
by analytic continuation: Gr(u>, k), seen as a function of u, can be analytically continued to 
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the upper half-plane, and for u n > one has 



G R (iu n , k) = -G E {u n , k) . 



(80) 



For simplicity, let's first examine the case where O = is dual to a massive bulk scalar 
field in the asymptotically- AdS background 



ds 2 = g^dx^dx" = e 2A(r) [-h{r)dt 2 + dx 2 } + e 2B(r) 



dr 2 
h(r) 



with the conformal boundary at r = and a black hole horizon with Hawking temperature 
T at r = rjj. The quadratic action for is 



S 



1 



d A x dr \J —g 



16ttG 5 

and the equation of motion following from it is 



g^d^d v <j> - > 2 2 



1 

—i 

2 



(82) 



□ - m 2 = . 



(83) 



Defining the Fourier transform of by 



u,k,r)= d A xe- l{t ^ t] (t){t,x,r) 



54) 



one can write 



as 



h! 



( AA' — B' + —\ d s 



e 2(B-A) 
h 



~— + k 2 + e 2A m 2 
h 



(w,k,r)=0. (85) 



For non-zero u, there are two linearly independent solutions to this equation: one for which 
m \th — r|~ Ju; / 47rT at r m rn, corresponding to 0-quanta falling into the black hole horizon, 
and one for which |r# — r| JaJ / 47rT , corresponding to 0-quanta coming out of the black 
hole horizon. The expectation of [20] was that the infalling solution should be the one 
related to the retarded correlator. At the conformal boundary, a generic solution behaves as 
m e ( - A_4 ' )A , where A is the more positive root of the equation 



A(A - 4) = m 2 L 2 



To go to Euclidean signature, one merely needs to replace —h(r)dt by h(r)dr in ([HI 
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and let r run from to (3. The Euclidean metric is 



ds% = g E ^udx^dx v = e 2A{r) [h(r)dr 2 + dx 2 ] + e 2B(r) 



dr 2 
h{r) 



(87) 



and the Euclidean action is 



S f 



1 



d A x dr J~g^ 
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If we define the Fourier modes of 0b by 



1 r 13 

ijsK, k ^) = pj dr I d 3 x e-^-^Mr, S, r) 



(89) 



then the equation of motion that follows from (IHHi) is 



u>\ p 2{B-A) / ,2 



E K,£,r) =0. (90) 



This equation is identical to ( |85l) except for the replacement uj 2 — >• — o;^. For non-zero w n , 
the near-horizon behavior of the two linearly independent solutions is <\>e ~ \th ~ r|' w "'/ 47rT , 
corresponding to regular solutions, and <pE ~ \th — r| _ '"™'' 47rT , corresponding to divergent 
solutions. The behavior of 0^ at the conformal boundary is the same as in the Minkowski 
case: generically, <pE ~ e^ A_4 ^, with A again defined as the more positive root of equation 



The Euclidean correlator (1791) can be computed from 



con-shell r jl I 
[<PE\ 



(91) 



where (j)o(u n ,k) = lim r ^ 0e(^n, fc, r)e^ 4_A ^ A and it is assumed that is regular at the 
horizon [131 [T5] . The on-shell action S E n ~ sh - cll [<j)E] can be obtained from integrating the first 
term in flHHj) by parts and using the equation of motion for <pE to set the bulk term to zero 
and keep only the boundary term. In Fourier space, the on-shell action can be expressed as 



con— shell P 



d 3 k 



J, 



E \VE\ 



(u n ,k,r) 



r=r H 



r=Q 



(92) 



Je [</>e] k, r) = (j) E (-u n , -k, r)he <9 r </> E (u; n , k, r) . 



We now introduce the notation for a solution of fl90|) that is regular at the horizon and 
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is normalized by 0o = 1- From (!9T|) it is straightforward to show that 

-. r=rn 

G E (u n , k) = (Je We S ] K, fc, r) + J E [4>e S ] r)) . (93) 

Note that the equation of motion ( 1901 for <pE as well as the boundary conditions are invariant 
under oj n — > — u; n , which implies that 

J £ [<^ cg ] (u, n , k, r) = J E [<^ cg ] (-w B , -k, r) . (94) 

So the four terms in fl93l) are equal in pairs. A further simplification occurs for u n ^ 0: 
starting from the fact that ~ |r# — r |l<^l/ 47rT at r « r# and using the definition of in 
( 1921) . one can show that [0^ g ] (o; n , fc, r) « \m — r\ 2 ^ n ^ 4nT — > as r — > r#. So there is no 



contribution from the horizon in ( 193]) . Consequently, for u; n 7^ 

1 



G E (u n , k) = — "— lim J E [<j)™ e } (u n , k, r) . (95) 

107TCt5 r^O 



In analogy to fl92|) . one can define 



1 /" dud?k 



r=r H 
r=Q 

AA—B ■ 



(96) 



J [4>] (u, k, r) = —<j)(—uj, —k, r)he d r <j)(u, k, r , , 



which can be obtained from fl82l) after integrating by parts and using the equation of motion 
( 183]) to set the bulk term to zero. We now introduce the notation <f) m for a solution of 
(|85p that is infalling and is normalized by setting lim r ^ (j) m (uj,k,r)e^~ A ^ A = 1 and <p out 
for a solution that is outgoing and has the same boundary normalization. From comparing 
equations (185]) and (190]) one can easily see that if u n > then ^ g (w n ) = (fi m (iu n ), while if 
u n < then (p r ^ e (u n ) = </> out {iuj n ) . Because the behavior near the boundary of 0^ g , m , and 
out is, by definition, the same regardless of the value of w or w„, it follows that 

iim Je ran K , £, r) = ( - •; [f ;] <**• § r » » •*><; w 

- lim r ^o J [0 ] (*u; n , fc, r) if u n < . 



Because ( 1971) and (195]) imply (180]) . it follows by analytic continuation^ that Gr should be 



5 The analytic continuation is unique if one further assumes that Gr(u>, k) decays fast enough at large 
complex u [31] . 
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given by 

Gr(0>, k) = -V lim n J l^] ( W > £ r ) • ( 9§ ) 
107TG5 r^O 

This prescription is exactly the one proposed in [2"U]FI 

The argument leading to fl98|) was based on the assumption of non- vanishing u> and uj n . 
For <jj = u n = 0, the equation of motion for <fi coincides with the one for One of 
its linearly independent solutions approaches a constant at the horizon, and the other one 
diverges logarithmically as a function of r# — r. In this case, the Euclidean correlator may 
pick up a contribution from the horizon, which would mean an added contribution to (1951) . 
But by continuity, we expect (|98|) to apply unaltered at u = if m is replaced by the 
solution which is regular at the horizon. 

The argument given above can be generalized to operators with higher spin, or to opera- 



tors that mix with one another. In the case of the shear viscosity computation in section l275j 
one can set O = T\2 and (f) = H\2- In the case of the bulk viscosity computation, there 
are some additional subtleties. To avoid complications present in the general case, in the 
remainder of this appendix we will consider only perturbations that do not depend on the x 
coordinates, which is to say we set k = 0. The fields Hn = H22 = whose dual operator 
O = jTi 1 we are interested in, couple to other fields, and so we must replace (f>(t, x, r) in the 
above discussion by h(t,r) = ^i/ 00 (t, r) Hn(t,r) H 55 (t, r)j, as in section 12721 Let's work 
out how the Green's function prescription (1981) generalizes to this case. 

To compute two-point correlators, we need only consider the quadratic part of the action 
for h. In terms of the Fourier transform 

h(u>,r)= r%^h(t,r), (99) 

the quadratic action is 

S = -J— / ^dr (d r h* T md r h + h* T h* T d r h + d r h* T hh + h* T kh) , (100) 

327TG5 J Z7T V / 

where m, b and k are the matrix- valued functions of uj and r given in (j22|) . In order to obtain 
an explicit final result, one would need to evaluate the contribution of the symmetric matrix 
G described below (I20dl) . We will not do this here, but we believe it is straightforward in 
principle. 



6 Note that if we wish to compute the advanced correlator we need only replace (j> m by <fi out . If we denote 
the result of following this alternate prescription by Ga(<jJ, k), then it follows from (|57)) that G J 4(*tJ n ,fc) = 
—Ge^u, k) for to n < 0, which is the analog of ([50)1 obeyed by the advanced correlator. 
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Integrating by parts and using the equations of motion that follow from fllOOp . one can 
see that (jHSJ) generalizes to 



J 



S 



h 



on— shell 



1 













h 


(u,r) 




l-oo 2tt 






r=0 
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(u, r) = h T (—uj, r) (md r + b) h(u, r) . 



(101) 



The Euclidean action Se is obtained from the real space counterpart of fllOOl) by performing 
the Wick rotation t — > —it and adding an overall minus sign. In Fourier space, the Wick 
rotation can be implemented by sending u to iuj n and replacing the integral over u by a 
discrete sum over n. One can thus write 



h, 



oon— shell 







Jk h 



[0J n ,r 



r=ru 
r=0 



32ttG 5 ^ 

n 

u n , r) (m E d r + b E ) h E (uj n , r) 



(102) 



where m.E{ui n ,r) = —m(iu) n ,r), and and kg obey similar relations. 

The solutions of the equations of motion that follow from (110 II) behave like infalling or 
outgoing waves near the horizon, namely h ~ \rjj — 7~| ±ia V 47rT f or r ~ Th- Here, in a slight 
abuse of notation, this asymptotic form is understood to hold for every component of h. It 
follows from this that at r pa rjj, the equations of motion for He have divergent solutions 
that behave as \th — r|~l w "l/ 47rT and regular solutions that behave as \th — r|' Wn '/ 47rT . 

To compute correlators of the operator dual to Hu, we must demand that Hqo — >■ 0, 
Hii — > 1 and H 55 —> as r — > 0. This normalization is not affected by Wick rotation 
and will be the same for the Euclidean fields He- It is useful to introduce the notation h^ e 
for a solution of the equations of motion for h E that is regular at the horizon and has this 
normalization near the boundary. We also will use h m and h ont to denote solutions of the 
equations of motion for h that are infalling or outgoing, respectively, and that are correctly 
normalized near the boundary. 

Equation (|93|) can be generalized to 



G E (uj n , k 



0) 



1 



32ttG, 



(uj n , r) + J E 



r=r H 



r=0 



(103) 



It is easy to see that J E 



(-u n ,r) = J E 



(u> n , r) and that J E 



h reg 



(uj n ,r H ) = 0, and 
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so 



G E (uj n , k = 0) 



1 



lim Je 



ireg 



(104) 



Finally, it is clear that h m (icu n ,r) = h T ^ g (u n ,r) for u n > and h out (iu n , r) = h r ^ g (u n ) for 
oj n < 0. The same reasoning as in the case of the bulk scalar field then shows that (|98|) can 
be generalized to 



G R (u, k = 0) 



1 



lim J 



(w,r) . 



(105) 
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The imaginary part of the retarded correlator can be expressed in terms of a conserved 
flux T . Indeed, (11001) is invariant under h — >• e lOL h and the conserved quantity associated to 
this symmetry is exactly 



h 



— Im J 



h 



(106) 



Since T is independent of r, it follows from (|98|) that ImGfl 



16vrG 5 ' 
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